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Theoretical approaches to particle propagation and acceleration in 
turbulent intergalactic medium 
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Abstract. Intercluster medium is expected to be turbulent with turbulence being superAlfvenic at large scales. Magnetic 
fields substantially modify the turbulent cascade when the turbulence reaches the scales at which the fluctuation velocity gets 
less than the Alfven velocity. At those scales it is possible to consider three cascades, of fast, slow and Alfven modes with 
little energy exchange between them. As Alfvenic and slow modes are anisotropic they marginally scatter and accelerate 
cosmic rays, while fast modes dominate the processes. However, in the presence of cosmic rays the turbulence is modified as 
cosmic rays transfer the energy of compressible motions (i.e. slow and fast modes) from large scales to the scale of cosmic ray 
Larmor radius. This results in generation of a new small-scale Alfvenic component which is not a part of the ordinary MHD 
cascade. This component does scatter and accelerate cosmic rays. In addition, magnetic reconnection in turbulent medium 
accelerates cosmic rays. The complexity of the interacluster turbulence calls for observational studies. A new technique 
Velocity Coordinate Spectrum (VCS) is particularly promising for studies of velocity fluctuations with a new generation of 
X-ray observatories. 
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^ 1. Do we expect IGM to be turbulent? 



A fluid of viscosity ly becomes turbulent when the rate of 
viscous dissipation, which is ^ v/L^ at the energy injec- 
tion scale L, is much smaller than the energy transfer rate ~ 
Vl /L, where Vl is the velocity dispersion at the scale L. The 
ratio of the two rates is the Reynolds number Re — VlL/v. 
In general, when Re is larger than Recrit ~ 10—100 the sys- 
tem becomes turbulent. Chaotic structures develop gradually 
as Re increases, and those with Re ^ lO'^ are appreciably 
less chaotic than those with Re ~ 10®. 

It is widely accepted that medium in clusters should be 
turbulent. Mergers of galactic subclusters may be one of 
the major energy injection mechanism (see Sarazin 2002, 
Brunetti 2003 and references therein). The details of the in- 
jection and the energy transfer are poorly understood. A crude 
picture includes the energy injection scale of 100-500 kpc and 
the injection velocity of the order of 10"^ km/s. 

A difficulty that one faces trying to understand turbulence 
is clusters is that the diffusivity is very different along and 
perpendicular to magnetic field. This is related to the fact that 
the mean free path of ions is substantial. For non-magnetized 
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intracluster gas the Re number is only marginally larger than 
Recrit and therefore formally is only just sufficient for initi- 
ating the turbulence. At the same time, even small magneti- 
zation makes the Reynolds number for the motions perpen- 
dicular to magnetic field very large (e.g. more than 10^°). 
This poses the question of what a reasonable choice of the 
Reynolds number is for intracluster medium. 

Strangely enough we do not have a good answer to this 
basic question. There are several processes that definitely in- 
fluence the effective diffusivity of intracluster plasma. First 
of all, compressions of collisionless plasma should result in 
various instabilities that would decrease the mean free paths 
of protons. In addition, bending of magnetic field lines does 
not allow mean free path of protons to be larger than the scale 
at which the turbulent velocity is equal to the Alfven velocity 
(henceforth, the Ia scale). Therefore the effective Reynolds 
number at the injection scale will be ^ 10'^ for the field of 1 
/zG field, which is the field that fills according to Enslin et al. 
(2005) 90% of the intracluster volume. 

Initially the turbulence is superAlfvenic and hydrody- 
namic motions easily bend magnetic field lines. Such turbu- 
lence is analogous to hydrodynamic one until the scales of the 
order of Ia = L{Va/Vl)'^ are reached. Assuming Vl — 10'^ 
km/s and L = 100 kpc we get 30 pc, if we adopt, following 



3 HOW DOES TURBULENCE SCATTER AND ACCELERATE COSMIC RAYS? 



Ensslin et al (2005) that 1 /xG field fills 90% of the ingra- 
cluster volume. The turbulence gets mangetohydrodynamic 
(MHD) for scales Z < /a as magnetic fields control fluid mo- 
tions. However, in both "hydrodynamic" and MHD regimes 
plasma effects continue to be important as they drain and re- 
distribute energy of compressible motions (see Schekochihin 
et al. 2005). 

2. What do we know about MHD turbulence? 

There have long been understanding that the MHD turbulence 
is anisotropic (e.g. ShebaUn et al. 1983). Substantial progress 
has been achieved by Goldreich & Sridhar (1995; hereafter 
GS95), who made an ingenious prediction regarding relative 
motions parallel and perpendicular to magnetic field B for 
incompressible MHD turbulence. An important observation 
that leads to understanding of the GS95 scaling is that mag- 
netic field cannot prevent mixing motions of magnetic field 
Unes if the motions are perpendicular to the magnetic field. 
Those motions will cause, however, waves that wiU propagate 
along magnetic field Unes. If that is the case, the time scale of 
the wave-like motions along the field, i.e. ~ 1\\/Va, is the 
characteristic size of the perturbation along the magnetic field 
and Va = B/ y^Anp is the local Alfven speed) will be equal to 
the hydrodynamic time-scale, l±/vi, where l± is the charac- 
teristic size of the perturbation perpendicular to the magnetic 
field. The mixing motions are hydrodynamic-like. They obey 

1 /3 

Kolmogorov scaling, vi a l_l , because incompressible tur- 
bulence is assumed. Combining the two relations above we 
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can get the GS95 anisotropy, ex Ij^ (or A:|| ex k_l in 
terms of wave-numbers). If we interpret l\\ as the eddy size in 
the direction of the local magnetic field, and lj_ as that in the 
perpendicular directions, the relation implies that smaller ed- 
dies are more elongated. The latter is natural as it the energy 
in hydrodynamic motions decreases with the decrease of the 
scale. As the result it gets more and more difficult for feeble 
hydrodynamic motions to bend magnetic field lines. 

GS95 predictions have been confirmed numerically (Cho 
& Vishniac 2000; Maron & Goldreich 2001; Cho, Lazar- 
ian & Vishniac 2002, 2003); they are in good agreement 
with observed and inferred astrophysical spectra (see Cho 
& Lazarian 2005). What happens in a compressible MHD? 
Does any part of GS95 model survives? The answer de- 
pends on the mode coupling. According to closure calcula- 
tions (Bertoglio, Bataille, & Marion 2001; see also Zank & 
Matthaeus 1993), the energy in compressible modes in hy- 
drodynamic turbulence scales as ^ Af^ if Mg < 1. Cho & 
Lazarian (henceforth CL03) conjectured that this relation can 
be extended to MHD turbulence if, instead of M^, we use ~ 
{SV)\/{a^ + Vl). (Hereinafter, we define Va = Bo/\/4^, 
where Bq is the mean magnetic field strength.) However, 
since the Alfven modes are anisotropic, this formula may re- 
quire an additional factor. The compressible modes are gener- 
ated inside the so-called Goldreich-Sridhar cone, which takes 
up ~ {SV^AjyA of the wave vector space. The ratio of com- 
pressible to Alfvenic energy inside this cone is the ratio given 
above. If the generated fast modes become isotropic (see be- 
low), the diffusion or, "isotropization" of the fast wave energy 



in the wave vector space increase their energy by a factor of 
^ Va/ {W)a - This results in 

6Ecomp SVaVa 

which suggests that the drain of energy from Alfvenic modes 
is marginal along the cascade. 

Our considerations above about the mode coupling can 
guide us in the discussion below. Indeed, if Alfven cascade 
evolves on its own, it is natural to assume that slow modes 
exhibit the GS95 scaling. Indeed, slow modes in gas pressure 
dominated environment (high /3 plasmas) are similar to the 
pseudo- Alfven modes in incompressible regime (see GS95; 
Lithwick & Goldreich 2001). The latter modes do follow the 
GS95 scahng. In magnetic pressure dominated environments 
or low /3 plasmas, slow modes are density perturbations prop- 
agating with the sound speed a parallel to the mean magnetic 
field. Those perturbations are essentially static for a <C Va. 
Therefore Alfvenic turbulence is expected to mix density per- 
turbations as if they were passive scalar. This also induces the 
GS95 spectrum. 

The fast waves in low P regime propagate at Va irrespec- 
tively of the magnetic field direction. In high /3 regime, the 
properties of fast modes are similar, but the propagation speed 
is the sound speed a. Thus the mixing motions induced by 
Alfven waves should marginally affect the fast wave cascade. 
It is expected to be analogous to the acoustic wave cascade 
and hence be isotropic. 

3. How does turbulence scatter and accelerate 
cosmic rays? 

The propagation of cosmic rays (CRs) is affected by their in- 
teraction with magnetic field. This field is turbulent and there- 
fore, the resonant interaction of cosmic rays with MHD tur- 
bulence has been discussed by many authors as the principal 
mechanism to scatter and isotropize cosmic rays. Although 
cosmic ray diffusion can happen while cosmic rays follow 
wandering magnetic fields (Jokipii 1966), the acceleration of 
cosmic rays requires efficient scattering. For instance, scatter- 
ing of cosmic rays back into the shock is a vital component 
of the first order Fermi acceleration. 

While most investigations are restricted to Alfven modes 
propagating along an external magnetic field (the so-called 
slab model of Alfvenic turbulence) (see Schlickeiser 2002), 
obliquely propagating MHD modes have been included in 
Fisk et al. (1974) and later studies (see Pryadko & Petrosian 
1999). A more complex models were obtained by combin- 
ing the results of the Reduced MHD with parallel slab-like 
modes . Models that better correspond to the current under- 
standing of MHD turbulence (see above) have been consid- 
ered lately. Chandran (2000, henceforth COO) considered res- 
onant scattering and acceleration by incompressible MHD 
turbulence. Resonant scattering and acceleration by the com- 
pressible MHD turbulence was considered in Yan & Lazarian 
(2002, 2004, henceforth YL02, YL04). There the following 
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result was obtained for the diffusion coefficients that gov- 
erned by Alfvenic scattering 
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where T[a, z] is the incomplete gamma function. This result 
was obtained using a tensor of magnetic fluctuations that was 
obtained in the CLV02 study. It provides scattering and ac- 
celeration rates orders of magnitude larger than those in COO 
for the most of energies considered. However, if energy is in- 
jected through random driving at large scale L ::g> k'^l^^ the 
scattering frequency, 

u = 2D^^/{1 - fx^), (3) 

are still much smaller than the estimates for isotropic and slab 
model. What does scatter cosmic rays? Our work in YL02 
identified fast modes as the principal agent responsible for 
CR scattering and acceleration. A study in YL04 showed that 
this is true in spite of the fact that fast modes are subjected to 
much more dissipation compared to the Alfven modes. In the 
next section we shall discussed another possibility proposed 
in Lazarian & Beresnyak (2006), namely, instabihties related 
to cosmic rays that inject Alfven modes at resonant scales 
(see §4). This possibility may provide a physical justification 
for the earlier calculations in Brunetti et al. (2004). 

The work above was done in relation to galactic cosmic 
rays. However, the problem of MHD turbulence with CR is a 
general one. Therefore in a recent paper Cassano & Brunetti 

(2005) considered fast modes as the principal component of 
CR acceleration in galaxy clusters (see also Brunetti 2004). 
A further work in this direction is in Brunetti & Lazarian 

(2006) . 

Resonant acceleration (we include Transit-Time- 
Damping acceleration to this category) is not the only 
process by which magnetized turbulence accelerates CR. 
For instance, the acceleration of cosmic rays by the large 
scale compressible motions was described in the literature 
rather long time ago (see Ruskin 1988). In this regime slow 
and fast diffusion hmit exist. The slow limit corresponds 
to the rate of particle diffusion out of compressible eddies, 
which is slower than the evolution rate of the eddies. On the 
contrary, in the fast diffusion limit particles leave the eddies 
before they turnover. Large scale compressions associated 
with anisotropic slow modes were used to accelerate Solar 
Flare CR in Chandran (2003) who considered slow diffusion 
regime. Fast diffusion regime was used in Chandran & 
Maron (2004ab) in application to CR acceleration in galaxy 
clusters. A comprehensive study of these processes has been 
done in Cho & Lazarian (2006). We found, first of all, in 
slow diffusion Umit the resonance scattering, that is a part 
and parcel of the slow diffusion process, is the dominant 
acceleration process. Then, we identified fast modes as the 
principal cause of non-resonant acceleration. In addition, we 
showed that weak turbulence (Galtier et al. 2005) may be 
important for cosmic ray acceleration in the fast diffusion 
limit. 



Alfvenic slab motions 




2x10 



scale, AU 



Fig. 1. Energy density of compressive modes and Alfvenic 
slab-type waves, induced by CRs, in galaxy clusters. The en- 
ergy is transferred from the mean free path scale to the CR 
Larmor radius scale. If mean free path falls below compres- 
sive motions cutoff, spectrum of slab waves becomes steeper 
(from Lazarian & Beresnyak, 2006). 



While the importance of fast modes for the scattering and 
the acceleration of CR is unquestionable, plasma and CR in- 
stabilities may make the actual turbulence in galaxy clusters 
more involved compared to a simple picture presented above. 
For instance, below we consider a robust CR instability that 
should modify the properties of MHD turbulence. 

4. How can cosmic rays modify turbulence? 

Studies of the backreaction of the energetic particles on the 
turbulence are usually limited by the damping of turbulence 
on energetic particles (see Brunetti & Blasi 2005, Ptuskin et 
al. 2005, Petrosian, Yan & Lazarian 2006). However, this is 
not the only effect of CR. For instance, Lazarian & Beres- 
nyak (2006) have considered a transfer of turbulent energy 
to small (but not dissipation!) scales that is mediated by CR. 
As CR present a coUisionless fluid, compressions CR through 
the compression of magnetic field preserve the adiabatic in- 
variant p\/B, where p_L is the CR momentum perpendicular 
to magnetic field. CR with the anisotropic distribution of mo- 
menta, i.e. with nonzero A = (px — P\\)Ip\\ subjected to 
an instabiUty, which growth rate can be estimated as 



Ti{k\\) = ujpi -AQ, 



(4) 



where il = eB/mc is a cyclotron frequency, m - proton 
mass, n is the density of plasma, iVpi is the ion plasma fre- 
quency, refering to n and ncnip > mfl/k) is the number 
density of cosmic rays with momentum larger than minimal 
resonant momentum for a wavevector value of k. Q is a di- 
mensionless numerical factor, depending only on cosmic ray 
power-law index a. 

The outcome of this instability is the direct transfer of the 
energy from mean free path of a CR, which is determined 
by magnetic scattering of the CR, to the CR gyroradius. The 
magnetic perturbations at the gyroradius scale are the most 
efficient in CR scattering and this decreases the mean free 
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path A and therefore A ^ vx/Va- At large scales the instabil- 
ity is limited by damping arising from the ambient Alfvenic 
turbulence (rp^crit in Fig. 1), while at small scales steep- 
ening limits the amplitude the perturbations. All in all, the 
interaction of CR with MHD turbulence results in an addi- 
tional slab-like small-scale component of Alfvenic perturba- 
tions that interacts with CR much more efficiently than the 
Alfvenic mode being excited by the energy injection scale. 

5. Can magnetic reconnection accelerate CR? 

Magnetic reconnection is a process that we expect to hap- 
pen routinely in magnetized turbulent fluid. This even more 
true for superAlfvenic turbulence when fluid motions bend 
magnetic fields easily. Indeed, magnetic reconnection should 
happen whenever magnetic field with non-parallel direction 
interact. 

Several schemes of magnetic reconnection are known 
(see Fig. 2). The most famous one is the Sweet-Parker 
(Parker 1957, Sweet 1958) recoimection, which naturally oc- 
curs when magnetic fields get into contact over long current 
sheets. This is the most natural arangement of magnetic field 
to get into within the turbulent fluid as random moving flux 
tubes collide and push their way through one another. The 
scheme provide ridiculously low rates of recoimection to be 
of astrophysical importance, however. Indeed, the recoimec- 
tion rate scales as VAi?m~^/^, where Rm = VAL/rimag 
is the Lunquist number, which for realistic values of mag- 
netic diffusivity rjmag is so humangous that any processing of 
magnetic energy via reconnection gets absolutely negligible 
for any astrophysical system (e.g. stars, interstellar medium) 
not to speak about intracluster medium with its much larger 
scales. 

The Petcheck (1964) scheme emerged as an answer to the 
evidence of fast reconnection that cannot be explained within 
the Sweet-Parker model. Within this scheme the reconnection 
is concentrated along thin filaments from which oppositely 
directed magnetic field lines diverge. It is clear at this moment 
that it cannot be sustained at large Rl for smooth resistivi- 
ties. Whether or not anomalous effects, e.g. those related to 
Hall term, can provide reconnection at rates comparable with 
the Alfven speed is hotly debated' (see Biskamp, Schwarz 
& Drake 1997, Shay et al. 1998, Shay & Drake 1998, Bhat- 
tacharjee. Ma & Wang 2001). We feel, however, that the issue 
of satisfying boundary outflow conditions is the most contro- 
versial element in applying Petscheck scheme to astrophys- 
ical conditions. If very special global geometry or magnetic 
fluxes, e.g. convex magnetized regions, is required to enable 
reconnection, then for a generic astrophysical case the recon- 
nection is slow (see Fig. 3). 

The turbulent reconnection model proposed in Lazarian 
& Vishniac (1999, henceforth LV99) deals with magnetic 

' For instance, the necessary condition for the anomalous effects 
to be important, e.g. for that the electron mean free path is less than 
the current sheet thickness (see Trintchouk et al. 2003) is difficult to 
satisfy for the ISM, where the Sweet-Parker current sheet thickness 
is typically much larger than the ion Larmor radius. 
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Fig. 2. Upper plot: Sweet-Parker model of reconnection. The 
outflow is limited by a thin slot A, which is determined by 
Ohmic diffusivity. The other scale is an astrophysical scale 
L » A. Middle plot: Turbulent reconnection model that ac- 
counts for the stochasticity of magnetic field lines. The out- 
flow is limited by the diffusion of magnetic field lines, which 
depends on field line stochasticity. Low plot: An individual 
small scale reconnection region. The reconnection over small 
patches of magnetic field determines the local reconnection 
rate. The global reconnection rate is substantially larger as 
many independent patches come together. 

field configurations with flat long current sheets. However, 
these sheets consist of small sheets related to the interaction 
of individual turbulent elements of magnetic flux. The out- 
flow within the model is limited not by the thickness of the 
current sheet, but by magnetic field wondering. As the result 
the rate of reconnection is proportional to the turbulence in- 
tensity. The model can explain flaring associated with recon- 
nection as well as high rates of reconnection that are required 
by both observations and theory (e.g. dynamo theory). Cal- 
culations in Lazarian, Vishniac & Cho (2004) confirm main 
points of the LV99 model (e.g. the rate of field line wonder- 
ing) and extend it to reconnection in partially ionized gas. 

There are several ways how magnetic reconnection can 
accelerate cosmic rays. It is well known that electric fields in 
the current sheet can do the job. For Sweet-Parker reconnec- 
tion this may be an important process in those exceptional in- 
stances when Sweet-Parker reconnection is fast in astrophysi- 
cal settings. For Petscheck reconnection only an insubstantial 
part of magnetic energy is being released within the recon- 
nection zone, while bulk of the energy is being released in 
shocks that support X-point. Therefore one would expect the 
shock acceleration of cosmic rays to accompany Petscheck 
reconnection. 

Similarly to the Petscheck scheme the turbulent reconnec- 
tion process assumes that only smaU segments of magnetic 
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Modification of Petscheck Reconnection 
by External Forcing 

Fig. 3. Upper plot: Petscheck reconnection scheme has a 
magnetic diffusion region (rectangular area at the tip of the 
magnetic field line bending) for which both the longitudinal 
and transversal dimensions are determined by the Ohmic dif- 
fusivity. To enable Vr ^ Va the model requires field line 
opening over the whole astrophysical scale involved, which 
is difficult to satisfy in practice. Lower plot: External forc- 
ing, e.g. the forcing present in the ISM, is likely to close the 
opening required by the Petscheck model. In this case the 
global outflow constraint is not satisfied and the resulting re- 
connection speed is V/^ <C Va- As the result the Petscheck re- 
connection cannot operate steadily and therefore cannot deal 
with the amount of flux that, for instance, an astrophysical 
dynamo would require to reconnect. 



field lines enter the reconnection zone and are subjected to 
Ohmic annihilation. Thus only small fraction of magnetic en- 
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ergy, proportional to (LV99), is released in the current 

sheets. The rest of the energy is released in the form of non- 
Unear Alfven waves that are generated as reconnected mag- 
netic field lines straighten up. Such waves are likely to cause 
second order Fermi acceleration. This idea was briefly dis- 
cussed in Lazarian et al. (2001) in relation to particle accel- 
eration during the gamma-ray burst events. In addition, large 
ampUtude Alfvenic motions in low /3, i.e. magnetically dom- 
inated, plasmas are Ukely to induce shocks (see Beresnyak, 
Lazarian & Cho 2005), which can also cause particle accel- 
eration. 

However, the most interesting process is the first-order 
Fermi acceleration that is intrinsic to the turbulent reconnec- 
tion. To understand it consider a particle entrained on a re- 
connected magnetic field line (see Fig. 4). This particle may 
bounce back and force between magnetic mirrors formed 
by oppositely directed magnetic fluxes moving towards each 
other with the velocity Vr. Each of such bouncing will in- 
crease the energy of a particle in a way consistent with the 
requirements of the first-order Fermi process. The interesting 
property of this mechanism that potentially can be used to 
test observationally the idea is that the resulting spectrum is 
different from those arising from shocks. Gouveia Dal Pino 
& Lazarian (2003) used particle acceleration within turbu- 
lent reconnection regions to explain the synchrotron power- 
law spectrum arising from the flares of the microquasar GRS 
1915+105. Note, that the mechanism acts in the Sweet-Parker 
scheme as well as in the scheme of turbulent reconnection. 
However, in the former the rates of reconnection and there- 
fore the efficiency of acceleration are marginal in most cases. 
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Fig. 4. Cosmic rays spiral about a reconnected magnetic field 
line and bounce back at points A and B. The reconnected re- 
gions move towards each other with the reconnection veloc- 
ity Vr. The advection of cosmic rays entrained on magnetic 
field lines happens at the outflow velocity, which is in most 
cases of the order of Va- Bouncing at points A and B hap- 
pens because either of streaming instability or turbulence in 
the reconnection region. 



6. How can we test the turbulence model? 

There are several ways that information about turbulence in 
intracluster medium can be obtained from observations. For 
instance, magnetic power spectrum has been obtained using 
Faraday rotation (see Ensshn 2004, Ensshn, Vogt & Pfromer 
2005). The issues of detectibility of velocity turbulence have 
been discussed in Sunyaev, Norman & Bryan (2003). Recent 
advances in the techniques of studies of turbulence via ve- 
locity fluctuations (see review by Lazarian 2004 and refer- 
ences therein) enable us to get spectra of turbulent fluctua- 
tions. The techniques described there can potentiaUy separate 
compressible and incompressible motions. For instance, a Ve- 
locity Coordinate Spectrum (VCS) technique can be used to 
extract turbulence statistics from the Doppler- shifted spectra 
even when the measurements are made over a Umited number 
of directions (even one) or the object is not resolved (Lazar- 
ian & Pogosyan 2006). Such measurements would already be 
possible, if not for the failure of the high velocity resolution 
instrument on board of the ASTRO-E2. An example of such 
a study that can be available with the forthcoming X-ray tele- 
scope ConsteUation X is presented in Fig. 5. 
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Fig. 5. An example of the velocity line profile (upper plot) 
and velocity spectrum (lower plot) for underlying Kol- 
mogorov turbulence that can be produced with observations 
using Constellation X forthcoming X-ray facilty with 1 hour 
exposure (Chepurnov & Lazarian 2006). 

7. What are the prospects of the field? 

The detailed description of the CR acceleration and scatter- 
ing is a challenging problem that requires coordinated efforts 
of both theorists and observers. It is clear that at the mo- 
ment we do not have an adequate description of MHD turbu- 
lence in the coUisionless fluid and this is an impediment for 
the progress. To get such a description one needs to account 
for different instabilities that affect magnetized coUisionless 
intercluster plasmas. This, together with the change of ef- 
fective Re fluid number as magnetic field and instabilities 
evolve, make the description of turbulence quite challenging. 
Together with turbulence driving the above processes deter- 
mine the efficiency of the CR interaction with magnetized 
turbulence. We know by now that if the energy is injected 
at large scales, the interaction of CR of low energies with 
Alfvenic part of MHD cascade is suppressed. This, however, 
may not be true if some part of Alfven modes is generated by 
small scale plasma instabilities. In particular, some of these 
instabilities may be due to CR themselves. In addition, re- 
connection within chaotic fields may accelerate CR directly. 
Therefore, with more observational data, clusters of galaxies 
Nevertheless, with more observational data, clusters of galax- 



ies may serve as a good testing ground for studies of the fun- 
damental processes involving CR. 
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